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Applying the method from recently developed fluctuation theorems to the stochastic 
dynamics of single macromolecules in ambient fluid at constant temperature, we es- 
tablish two Jarzynski-type equalities: (1) between the log-mean-exponential (LME) 
of the irreversible heat dissiption of a driven molecule in nonequilibrium steady-state 
(NESS) and lnP"'^'^''(a;), and (2) between the LME of the work done by the inter- 
nal force of the molecule and nonequilibrium chemical potential function ii^^^^{x) 
= U{x) + A;BTlnP""""(x), where P"^^"(x) is the NESS probability density in the 
phase space of the macromolecule and U (x) is its internal potential function. ^ = 
J jji^'^^^ (x) P"^^^^ (x)dx is shown to be a nonequilibrium generalization of the Helmholtz 
free energy and A^' = AU—TAS for nonequilibrium processes, where S = —ks J P{x) In P{x)dx 
is the Gibbs entropy associated with P{x). LME of heat dissipation generalizes the 
concept of entropy, and the equalities define thermodynamic potential functions for 
open systems far from equilibrium. 
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Introduction 

One of the important concepts in equilibrium thermodynamics is chemical potential. In an 
equilibrium, the chemical potentials of all participating chemical components are equal. Chem- 
ical potential also plays an important role in isothermal, driven, single molecular motor systems 
in nonequilibrium steady-state (NESS) lU. Potential functions for stochastic systems has been 
investigated previous by Graham and Haken |l2l|3l|4|| and recently by Ao and coworkers El ID. In 
all the above works, the logarithm of the probability distribution in a NESS that is analogous to 
Boltzmann's entropy has been identified as a nonequilibrium potential function (NePf). 

In classical mechanics and field theory, introducing a potential function is always associated 
with integration of a conservative force field. Thus, for the NePfs introduced above, one necessarily 
seeks an "integral interpretation". This is the objective of the present work which shows that, in a 
statistical sense, the integration is provided by the recent developed fluctuation theorems (FTs) and 
Jarzynski-type equalities in NESS EUHJ. FTs were first developed in the context of microcanonical 
chaotic dynamical systems ||9l[T0||. It has also been studied in parallel for stochastic canonical 
dynamics ifTTl fT2ll . In this paper, we first present our results via a simple example: a 3-state, 
driven cyclic reaction widely found in biochemistry. Then the more general derivations based on 
the polymer representation of a macromolecule are given. Useful mathematics can be found in the 
Appendices. 

NePf and log-mean-exponential of heat and work: an example 

Fig. 1 shows a driven cyclic reaction of a single enzyme in aqueous solution at constant tem- 
perature T, sustained by constant concentrations of cofactors D and E. This can be rigorously 
accomplished in a test tube with a feedback regenerating system, or approximately with a large 
excess of the cofactors [13|. In a living cell, such NESS reactions, e.g., kinase catalyzed phos- 
phorylations and GAP-activated GTP hydrolysis [14J, are ultimately maintained by the working of 
mitochondria which constantly regenerate ATP from ADP 

The chemical potential for a species X {= D or E) is fix = l^'x + ^bT In cx where cx is the 
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concentration for the species X. Then it is easy to show that 

Ahed = fJ'D - fJ'E = kBTln- — — . (1) 

The chemical potentials for the single macromolecule are nx = /^x + ^bT In px where px is the 
probability of the macromolecule being in state X (= A, B, or C). If the reaction is in a closed 
vessel, it eventually reaches an equilibrium and A^ed = ^ ^'^e/^d = ^~fcffc~4t~^ ^ c'^^ed/'^bT _ 
However, if one controls the reaction such that both concentrations of D and E are maintained 
at some constant values, say cd and ce, then Eq. [l]gives the amount of work done to the system in 
order to sustain the NESS. The NESS flux 

J = kipT" - k-ipT' = k2pT' - k-2Pc^"" = ksp'l?'' - k-3pT' 

_ kik^k^cp - k_ik%k_3CE 

kik2 + kik'i + k2k'i + kik_2 + ^2^-3 + ^3^1 + ^-1^-2 + ^-1^-3 + ^-2^-3 ' 



We thus have J x AfiED, the work done by the external agent per unit time, > and being zero if 
and only if when A/jed = J = 0- The analogy between this and power = current x voltage, from 
electrical circuits, should be noted. The work also is equal to the sum of the dissipated heat from 
each every isothermal reaction: A/jed = Qi + Q2 + (I3 (see Fig. 1). These results reflect both the 
First and Second Laws of Thermodynamics in a single macromolecule far-from-equilibrium [ 16|. 
By far from equilibrium, we mean in general for a reversible reaction X ^ y in a NESS network, 

k- 

J = k+cx — k-Cy = J+ — J- and Afiyx = ksT \n{k+cx / k^cy) = kBThi{J+/ J_) are not linearly 
related ifTT^. However, if J ^ J+, then one has A/i oc J. 

So far, all the materials are from standard chemical thermodynamics. (Though seldomly it has 
been presented in this way.) To introduce the new results from the recent work on FTs ifTTlfTllfTSll . 
we first define the irreversible heat dissipation following a stochastic trajectory in a NESS, Qirr{t)- 
Let the single macromolecule having the trajectory xq, xi, Xn (Xi = A, B, or C) with transition 
from state Xi^i to state Xi at ti (to = < ti < ... < tn < t), and defining 

n 

Qirr{t) = - ^ A(^^^_^^^ (3) 
i=l 
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where AfiAB = Qi, ^f^BC = 12, and AficA = 13, all positive. Qirrif) is called irreversible heat 
since the dissipation here are truly irreversible: Qirr{t) = in an equilibrium steady-state. There 
is, however, a fluctuating reversible heat exchange between the system and its environment even in 
equilibrium. The latter term is not contained in the Qirr{t)- In a NESS, for each every reaction A/x 
= Ah° — TAs° — T As concentration, hcucc the total heat dissipation = Afi — Ah° = —TAs, where 
entropy change As includes macromolecular structural entropy change As° and concentration 
entropy change As concentration- According to the recent work of U. Seifert ifTHll . A/i/T is the total 
(irreversible) entropy production. As is the entropy change, and Ah" /T is the exchanged heat 
which can be positive or negative. If the reaction is in an isothermal equilibrium, A/x = and 
fluctuation in heat exchange is equal to the enthalpy fluctuation along a stochastic trajectory. This 
exact compensation between fluctuating entropy and enthalpy was first discussed in [il9i . It is a 
direct consequence of microscopic energy conservation. 

Applying the methods developed for FTs ifT^lfTHlEUlETl . it can be shown that 

hipY' = In (5x,x„e-^-'-(*)/'^^^^ (4) 

where (• ■ ■ ) is ensemble average over all the possible trajectories with initial distribution p"'^^* for 
Xq, and 5jj is the Kronecker-delta function for discrete i and j. This is a new result. In the light of 
Jarzynski's equality and the work of Hummer and Szabo ||71|22|, Eq. Infirmly establishes Inp'^^^'^ 
as a NePf of the NESS. More specifically, these previous studies have shown that the integration 
of work (or heat), though itself is not path-independent and a function of the state, its LME is 
a function of state. Jarzynski in fact has shown that the LME of work is the free energy of the 
system. Now Eq. |4l shows that the LME of heat dissipation is a generalization of Boltzmann's 
entropy. Note that for equilibrium steady-state, Qirr{t) = and Eq.|4lbecomes a trivial equality. 

As the celebrated Jarzynski's equality, Eq. |4|is valid for all finite t. We note that Eq. lUdoes 
not provide a novel way to compute p"^'^'^ since computing heat dissipations, the g's, requires 
knowing p'^'^^^ a priori (Fig. 1). However, if we replace Qirr{t) by the internal work Wint{t) = 
- EILi = + /^xo' then we have 

fc^Tln (5x,.„e-^-W/^'^^)''° = (/i^ + kBT\npx{t)) - « + kBTlnp^M) , (5) 

where (■ • • )^'' denotes ensemble average over all possible trajectories which start at xq. The right- 
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hand-side is the chemical potential difference between X at time t and Xq at time zero. This result 
resemble very much the Jarzynski equality, in which the two end states are equilibrium states, 
while here the two end states are two configurations of a NESS. 

NePf based on stochastic macromolecule mechanics of driven single molecules 

The above results for the 3-state driven system can be generalized to any single macromolecules 
in aqueous solution with stochastic dynamics, either in terms of discrete master equations ll23ll or 
continuous Brownian dynamics ll24lfT9l . The discrete approach is a tradition of chemical kinet- 
ics. For a single macromolecule, the kinetics consists of unimolecular reactions with first-order 
rate equations and pseudo-first-order rate equations when cofactors are involved. The continuous 
Brownian (Langevin) approach to the dynamics of single macromolecules began with Kramer's 
theory ll25ll and has been extensively studied in polymer dynamics ll26l . We shall develop the 
general theory of NePf based on the Brownian dynamic formalism below. 

Based on the polymer dynamic approach, we have recently developed a nonequilibrium ther- 
modynamic theory for overdamped single macromolecules, such as a motor protein, in aqueous 
solution, which we called stochastic macromolecular mechanics lfT9l l27l l28ll . We have shown 
how thermodynamic concepts such as entropy production rate, heat dissipation rate, and Onsager's 
thermodynamic force are rigorously developed from the dynamic equation lfT9l . Furthermore, 
we have introduced mesoscopic, trajectory-based fluctuating energy Ut = U{Xt) and entropy 
T( = —ks In P{Xt) where Xt is a stochastic trajectory of a macromolecule in its 3A^-dimensional 
phase space and P{x) is the probability density. We observed that in an isothermal equilibrium 
Ut — TTt = along a stochastic trajectory (Recall the above example.) Beyond equilibrium 
steady-state, we have suggested that = U{x) — TT(x) as the nonequilibrium chemical po- 
tential function for the macromolecule in its driven NESS. Introducing the fluctuating entropy Tf 
in terms of the P{x), however, is not completely satisfying since it requires the probability density 
being solved from a Fokker- Planck equation a priori. Also from an experimental standpoint as well 
as a molecular simulation standpoint, this definition is not very useful since it can not be obtained 
along stochastic trajectories Xt. The fluctuating energy in Eq. |5l however, is computable from a 
stochastic trajectory (The statistical accuracy, of course, is an important but separated issue.) 
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In this section we establish the NePf, — InP(x), and the nonequilibrium chemical potential, 
IJ,{x), through FTs-like equalities (Eqs. |lland|51) in temis of stochastic dynamics of a single macro- 
molecule in aqueous solution ifT^ 

7]dXt = F{Xt)dt + y^2k^ dBt. (6) 

In Eq. |6lF(x) is decomposed into an internal potential part and an external driving part: F(x) 
= —VU{x) + G{x). ?7 is a frictional coefficient and we assume it is a scalar for simplicity. Bt is the 
standard Brownian motion which represents the collision between the atoms in the macromolecule 
and the solvent lfT9ll26l . The probability density function for Xt is the solution to a Fokker-Planck 
equation 

dP{x, t)/dt = -V ■ 3{x, t) = {kBT/7])V^P - (l/r/)V ■ (FP) . (7) 

All FTs focus on the work done by various stochastic forces. There are three forces in our 
problem: the total force F, the external force G, and the Onsager's thermodynamic force H = 
F -keTV In P{x,t) IH: 

Wtot{t)= [ F{X,)odXs, W,,t{t)= [ G(X,)orfX„ Q,rr{t)= [ n{Xs,s)odX,. (8) 
Jo Jo Jo 

These stochastic integrals are called total work, external work, and irreversible heat dissipation/total 

entropy production, respectively. We use 'o' for Stratonovich and '■' for Ito integrations. 

If the macromolecule is in a closed vessel, i.e., G(x) = 0, then both F and H are conservative 
forces with potential functions U (x) and U (x) — ksT In P{x, t) + Const., the internal energy and 
Helmholtz free energy (generalized), respectively ll29ll . More importantly, in equilibrium, i.e., the 
stationary state of the closed system, P^'^(x) oc Q-u{x)/kBT ^ n = as expected: The Onsager's 
thermodynamic force is zero. Then Wtot{t) is the fluctuating potential energy and Qirr{t) is the 
fluctuating free energy. Qirr{t) = when the macromolecule is in an equilibrium lfT9ll . 

It is easy to show that 

{Qirr){t) = j ^^n(a;) ■3{x,s) =r] ds j dxP~^3^{x, s), (9) 

is the mean entropy production which is always non-negative. 

{Wtot){t) = j^ds j dxF{x) ■ J(x, s), (10) 
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however, is the mean heat dissipation which can be negative in transient relaxation. Wtot{t) being 
the total heat dissipation has been discussed by several authors |I3ni[l7l[T2|. Finally, since 

We.t{t) = Wtotit) - Wintit) = Wtotit) + U{Xt) - U{Xo), (11) 

we see that the work done by the external force is equal to the heat dissipation plus the increase 
in internal potential energy. This is the trajectory-based energy conservation. In a NESS the 
internal potential energy fluctuations and entropy fluctuations are both stationary, i.e., E [U (Xt)] 
and E[\nP{Xt)] are constants, hence {Wext)it) = {Wtot)it) = {Qvn){t). The work done by 
the external force is completely transformed into dissipated heat, which precisely equals to the 
isothermal entropy production, in the NESS [31 1. 

In terms of irreversible heat Qirr{t), we now introduce 

t) = (^6{x - Xje-^-W/^'^^^ (12) 

where (■ ■ ■ ) denotes ensemble average over all possible Xt starts with R{x, 0) = P^^^^{x). It then 
can be shown mathematically that the following equality holds true for all time t (see Appendix 
A): 

lnP""^"(x) = In (5{x - Xt)e-'3-W/'=sTy (^3) 

which is our main result. It proves Eq. |H The trajectory-based log-mean-exponential (LME) of 
the irreversible heat dissipation is the NePf, as previously suggested by Graham (31, Ao [3|, and 
Qian [1|. Note that the "Boltzmann entropy" in Eq. \T^is defined through multiple trajectories 
with initial distribution p^^^'^. Alternatively with arbitrary initial distribution, Eq. [13] holds only 
for t — > oo. 

We now consider Wint{t), the work done by the internal potential force —VU. Wint is related 
to Wtot and and Wext (Eq. ^ and its ensemble average is zero in a NESS. The exponential average 

Q{x, t\y) = (^5ix - x,)e-"^-*W/^^^)' (14) 

where (■ ■ ■ )^ denotes ensemble average over all Xt which start at Xq = y. Noting that Wmtit) = 
U{Xq) — U (Xt), it is easy to see that 

ksT In Q(x, t\y) = U{x) + keT In P{x, t\y) - U{y) = t) + const. (15) 
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is, off by a constant, the nonequilibrium chemical potential at time t. In other words, the trajectory - 
based LME of the internal potential energy is the nonequilibrium chemical potential. When t 
oo, we have 

fcfiTln (5{x - Xi)e-^-*W/'=^^y ^ = U{x) + A;BTlnP"^^"(a;) - U{y). (16) 



Eqs. [T51 and [161 for transients and the NESS of open driven systems, our second main result, are 
parallel to Jarzynski's equalities ||71|22| for time-dependent but detail-balanced systems. 

The mathematics behind Eq. [151 is trivial since the Wmt is the work done by the internal 
potential force. In fact, to obtain P(x, t), the simplest way is by mathematical definition P(x, t\y) 
= {5{x — Xt)Y . We emphasize that the physically novel result is Eq.[T3l which defines the In p^'-^^^ 
as a thermodynamic potential. The physical meaning of Eqs. [151 and [161 then, is to establish 
equalities between other thermodynamic potentials and the LME of other types of work in an open 
system. In this spirit, we have a third relation of such (Appendix B): 

keT In (^5{x - Xt)e-^'^'^^^/^^^y = t) + const. (17) 

where the constant is independent of y and t; It is a function of x. 

NESS thermodynamics of single driven macromolecules: a theory 

With the chemical potential fx{x, t) defined, the free energy of a single macromolecule \l'(t) = 
J fi{x,t)P{x,t)dx has the fundamental property expected from a nonequilibrium generalization 
of Helmholtz free energy |[T9ll : 



d f 



= E[We,t{t)]-h,-Tj^S (18) 

where Cp = {Qirri't))/t = rj J ||JpP~^(ix is total entropy production rate, hd = {Wtotit))/t = 
J F ■ Jdx is the heat dissipation rate ifT^ . and S = —ksT / P In Pdx is entropy. The first term in 
([T51 is the work done by external work - energy input. The sum of the first two terms is the energy 
stored in the system - potential energy increase, AU. Hence, the relation A\l/ = AU — TAS hold 
true even for transient driven nonequilibrium processes. In a NESS, both d^/dt and dS/dt are 
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zero and the energy input are transformed entirely into dissipated heat. Conversely, can 
be written in terms of functional differentiation, as in classic thermodynamics: 

= — — . (19) 

t) also plays an important role in the transient dynamic relaxation to a NESS. One can intro- 
duce chemical potential difference t) — ^^^'^^{x) and its average 

H{t) = J Pix, t) t) - /i--(x)) dx = J P(x, t) In ( /il'^^]^ ) dx. (20) 

This is a Lyapunov function for the stochastic dynamics (Eq. iTj: H(t) > and {d/dt)H{t) < 
11241 1291 . Therefore, a transient process relaxing to a NESS minimizes, on average, the chemical 
potential difference with respect to its NESS. 

Hatano and Sasa have recently developed a theory for steady-state thermodynamics of 
stochastic systems. The relation between their work and ours ifTTl fT9ll is as follows. The key 
thermodynamic quantities in their theory are total heat dissipation Wtot{t), housekeeping heat dis- 
sipation, which is precisely our Qirr{t) in Eq. HI and mesoscopic entropy = —kB\nP^'^^^{x). 
As we have shown, the mean housekeeping heat (Qirr) is exactly the entropy production rate 
(Cp) in our work ifTTl fT9ll . this relation is further explored in ll32l . Their irreversible force b(x) 
= G{x) — V;u"^^'*(x) is exactly 11, the Onsager's thermodynamic force. We derived an entropy 
balance equation in a NESS dS/dt = Cp — E^y/totit)] which they called excess heat Qex- Their 
study, however, has gone beyond simple NESS systems and their main result is concerned with 
time-dependent driven system, i.e, F in Q is also an explicit function of t. In parallel with and 
treating T(x, t) as a time-dependent potential energy function, they introduced a charging process 

and found an interesting generalization of Jarzynski's result to driven systems. The mean (Y{t)) 
is a generalization of the relative entropy H{t) in (EUt . The present work, however, is limited to 
time-independent F(a;) and is concerned with establishing the concept of NePf and nonequilibrium 
chemical potential in terms of the LME of certain type of work in a nonequilibrium open system. 

Discussions 
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Stochastic NESS circuit theory. An analogy between the nonequilibrium steady- state thermo- 
dynamics and electrical network theory exists, fi represents potentials, V/i represents potential 
gradient, J represents flux, and V ■ J = is Kirchhoff 's current law. Tellegen's theorem then says 
that V/i, is perpendicular to J in a NESS: 



In the case of a system being driven at its boundary, F(x) = — Vf/ (x) forx&Q and the nonequi- 
librium driving force is only imposed on the boundary dVl: 



is analogous to Kirchhoff 's potential law. For x G dQ, J 7^ is called flux injection and fi{x) 7^ 
is called potential clamping |I3T1 ESI El • If there is neither flux injection nor potential clamping, 
then the system (x E ^l) is closed and it reaches equilibrium with J = Vjj. = 0. 

Graham-Haken decomposition. In the present work, the decomposition of nonconservative 
force F(x) = — Vf/(x) + G{x) is based on the physical reality of open and closed systems: Any 
open system has a unique corresponding closed system with internal potential U{x). Knowing 
this correspondence greatly enhances the thermodynamic understanding of the open system (see 
dEl). Mathematically, however, Helmholtz-Hodge decomposition F = — V0 + 7 is not unique. 
Graham and Haken [2| have shown that if one chooses (p = — ksT In P^'^'^^{x), then 7 satisfies 
V ■ 7 — V0 ■ 7 = 0. In terms of this particular decomposition, our in Eq. [151= const., and more 
importantly, Q^^{x,t) = J Q(x,t\y)dy = 1. That is, the exponential average of Graham-Haken 
"internal work", over an initially uniform distribution, is 1 for any finite time t. The significance 
of this result is that one, having experimental or simulated realizations of many trajectories of Xt 
with finite time duration (0 < t < t*), can devise an inverse algorithm to construct such 
that 



Finally, we notice that if one chooses the Graham-Haken decomposition, then by the result in 
footnote 3, we have 




(22) 




(23) 




(24) 



{6iy-Xt)e-^o^^^s)odx. 



(25) 
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in the limit of t ^ oo. 

Umbrella sampling in stochastic dynamics. Umbrella sampling is an important method in 
molecular simulations ll35l . Our result on nonequilibrium chemical potential also has an intimate 
relation to a more general mathematical result which has applications for umbrella sampling. Con- 
sider an arbitrary function V{x) and integral /y(t) = — /g W{Xs) o dXs, then 

{^{5{x - X,)e-^-«/'=^^)'e(^(J^)-^(^«/'=^^) = P^''\x) (26) 

in which the inner (■ ■ ■ )^ is an averaging over all Xt with Xq = y and the outer (■ • • ) is an 
averaging over all Xq with the stationary distribution of X. The result is the stationary distribution 
which can be obtained from sampling of Xt weighted toward those trajectories which move along 
with the gradient of V^(x). 

Conclusions. Open, driven chemical and molecular systems far from chemical equilibrium 
have been studied from several different points of view in statistical physics. In terms of stochastic 
processes, these investigations include "nonequilibrium potentials" C] |4l [Si, "dissipative struc- 
tures" EH, "cycle kinetics" EH |33 EH, "Brownian motors" (HIISI, "stochastic resonance" 
ll4m . and "fluctuation theorems" ifTTl fT2l [TSl fTSll . The result we obtain in this work, in terms 
of a Jarzynski-type equality, clearly demonstrates an intimate connection between the fluctuation 
theorems and the nonequilibrium potentials. Combining with previous work which has shown re- 
lations between Brownian motors and cycle kinetics PTIl . between cycle kinetics and fluctuation 
theorems ifTTI . between Brownian motors and excess entropy ll42l . and between stochastic reso- 
nance and cycle kinetics 143^, a consistent and comprehensive statistical thermodynamic theory for 
single macromolecules far from equilibrium emerges. Finally but not least, a variational principle 
for NESS in terms of LME of heat dissipation seems to be in sight. 

I thank P. Ao, Z.Q. Chen, R. Fox, G. GaUavotti, R. Graham, C. Jarzynski, S.-i. Sasa, J.M. 
Schurr, U. Seifert, Z.F. Shao, and X.S. Xie for helpful discussions. 
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Appendices 

In stochastic mathematics, Langevin equations such as in Eq.l^are known as stochastic differ- 
ential equations [44|. Because the random force dBt, a solution to the equation is stochastic and 
has to be characterized in terms of a probability density function P(x, t) which satisfies the corre- 
sponding Eq.|7] For arbitrary force K(x, t), let t) = — Xt)e^^o ^{Xt,t)odXt/kBT^ where 
(■ ■ ■ ) denotes ensemble average and let (j){t,x) = (^e~ -^o T<-{Xt,t)odXt/kBTY where (■ ■ • )'^' denotes 
ensemble average over all possible trajectories (solutions to Eq. ^ which start at x, the Cameron- 
Martin-Girsanov formula iH^ states that (j){t,x) satisfies the equation dcjy/dt = C*(j) and (f){x,t) 
satisfies the equation d(j)/dt = C(f), where the elliptic operator C = {kBT/ifjV'^ — {l/r])\F ■ V + 
V-F-2K-V-V-K - (1/A;bT)K ■ (K - F)], and is the adjoint of C. 

Now for arbitrary scale q{x,t), let ip{x,t) = — Xt)e"^''^^^'''*^'^'^/'^^^) and ip{t,x) = 
(^Q- lo q(Xs,s)ds/kBTy ^ the Feymann-Kac formula [44] states that^/'(t, x) satisfies the equation dip / dt 
= C*il) and ip{x,t) satisfies the equation dip/dt = Cip, where the operator C = {kBT/ifjV'^ — 
( l/r/) [F ■ V + V ■ F] — ( 1/ kBT)q{x, t) . These two adjoint relations are behind the work of Jarzyn- 
ski Q and that of Hummer and Szabo In those studies, the work W{t) is defined through an 
integral /J F{Xs, s)v{s)ds where F{x, s) is a time-dependent external force and v{s) is a deter- 
ministic rate. Hence Feymann-Kac, instead of Girsanov, formula was used. 

Appendix A. R{x,t) = {6{x - Xt)e"'3-' W/'=sr^ ^^t^ R{x,0) = P^'^'^x). That is the above 
K = n in Eq.lHl Hence R{x, t) satisfies the partial differential equation dR{x, t)/dt = CR{x, t) 
where 

£ = -^V^ - - [F ■ V + (V ■ F) - 2n ■ (V - VlnP"'^^*)] . (27) 

T] T] 
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It is easy to verify that P"'^'**(x) is the solution to this equation for all time t. Hence we arrive at 
Eq.ini 

Appendix B. If the F in Eq. |6lhas F = Fi + F2, then applying the above result we have 

{6{x - Xt)e-^o'P^i^s)odXs^y ^ (^^(^y _ x^)e-Io'P2i^^)''<i^sy_ ^ direct consequence of this result 

is a symmetry relation {6{x - Xt)e~^^'°*(*))^ = {5{y - Xt)e"(^~^)'^*°'(*))^. This is an alternative 
statement of the FTs. 

Now let Z{y, t) = {6{x — Xt)e~^'=^'(*)/''^^)^. As a function of y, it satisfies the Kolmogorov 
backward equation dZ{y, t) /dt = MZ{y, t) where the operator M = {kBT/r])V'^ - (V^) [F- V + 
(V ■ F) + 2{VU) ■ V + V^U] + {l/T]kBT)VU ■ G = e^/^'^^^oe^^/'^^^. £0 is the right-hand- side 
of Eq.|7l Because Wext + Wmt = Wtot, we arrive at Eq. [l7| 
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^MED=qi+<i2+q3 



Figure 1: A single macromolecule can be in states A, B, and C. The transition from i? to C is 
accompanied by the turnover of a cofactor between D and E via the second-order fcg and k°_2. 
We denote pseudo-first-order rate constants /ca^D = ^2 and k°_2'^E = ^-2- The total chemical 
driving force for the cyclic reaction is AfiED given in Eq. [U gi = ksT \n{pAki/pBk-i), q2 = 
ksT \n{pBk2/pck-2), and qi = kBT\\i{pckz/pAk~.z), are the heat dissipation associated with steps 
1, 2, and 3, respectively. The equality in the box reflects the conservation of energy. 
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